Critical Level-Spacing Distribution for General Boundary Conditions 
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It is believed that the semi-Poisson function P(S) = 4Sexp(— 2S) describes the normalized dis- 
tribution of the nearest level-spacings S for critical energy levels at the Anderson metal-insulator 
transition from quantum chaos to integrability, after an average over four obvious boundary con- 
ditions (BC) is taken (Braun et al 1]). In order to check whether the semi-Poisson is the correct 
universal distribution at criticality we numerically compute it by integrating over all possible bound- 
ary conditions. We find that although P(S) describes very well the main part of the obtained critical 
distribution small differences exist particularly in the large S tail. The simpler crossover between the 
integrable ballistic and localized limits is shown to be universally characterized by a Gaussian-like 
P(S) distribution instead. 
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Chaos in quantum systems is known to appear in the 
level statistics of the spectral fluctuations for stationary 
energy levels obtained from Hamiltonians defined in ap- 
propriate Hilbert spaces @- In the quantum chaotic 
regime the loss of most symmetries is manifested via 
the appearance of level-repulsion between nearest neigh- 
bor energy levels which implies a highly correlated spec- 
trum |3J. Integrable systems display, instead, uncorre- 
cted spectra with random levels which allow spectral de- 
generacies as unambiguous signatures of high symmetry. 
The solution of Schrodinger equation for billiards having 
zero potential inside and infinite outside a bounded do- 
main manifests integrable or chaotic quantum behavior 
depending on the shape of the domain boundary. This 
quantum crossover corresponds to the classical transition 
from integrability to chaos in billiards 4]. Wave chaotic 
behavior has recently been observed experimentally for 
optical billiards |5|. 

Disorder in quantum systems, e.g. a random poten- 
tial for electrons in a lattice, is known to cause an An- 
derson metal-insulator transition from extended to local- 
ized wave functions 0. This is a true quantum phase 
transition from chaos to integrability with the extended 
wave functions of weakly disordered metals correspond- 
ing to quantum chaos and the localized wave functions of 
strongly disordered insulators to integrability. By adding 
weak disorder to a perfect system a ballistic to chaotic 
quantum crossover also occurs, similarly to quantum bil- 
liards. Only for stronger disorder one can obtain the 
metal-insulator transition from diffusive extended to in- 
tegrable localized states. One dimension is exceptional in 
this respect since no intermediate quantum chaotic be- 
havior can be seen with a crosover straight from ballistic 
to localized behavior. 
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The critical region at the metal-insulator transition be- 
tween diffusion and localization has also been explored 
for level statistics in terms of the nearest neighbor level- 
spacing distribution P(S) 0>@]- The disordered problem 
has certain advantages, such as the existence of many 
energy levels at the transition and a statistical ensemble 
over disorder which complements the usual energy ensem- 
ble. The corresponding critical P(S) turns out to be a 
scale-invariant hybrid of chaotic Wigner-like linear small 
S and an integrable Poisson-like exponential for large 
S. The main surprise, however, in these studies was the 
generic sensitivity of the obtained critical level statistics 
to boundary conditions (BC) This is now well under- 
stood since it simply reflects the nature of critical wave 
functions, which are in between extended and localized, 
described by scale-invariant fractal distributions. Their 
filamented structure, closer to one-dimensional, brings 
little overlap at the boundaries so it is accompanied by 
extreme sensitivity in boundary conditions. This is the 
reason for the corresponding generic dependence of criti- 
cal level statistics to BC. In this Letter we generalize the 
argument of Ref. Q by including many possible BC via 
the introduction of all Aharonov-Bohm fluxes and per- 
form an integration over them. 

Our aim is to compare the obtained level-spacing dis- 
tribution at the metal-insulator transition with the sim- 
ple semi-Poisson form P(S) = 45exp(-2S") 9]. Al- 
though a similar hybrid distribution describes, at least 
partially, the level-spacings in various metal-insulator 
transitions in disordered systems including a system with 
spin-orbit coupling |10| , the precise semi-Poisson form of- 
ten fails to describe criticality and is not widely accepted 
as the universal critical distribution. We ask the ques- 
tion whether a semi-Poisson might be the appropriate 
critical distribution to describe averages over all possible 
BC. For this purpose we set up and diagonalize Hamil- 
tonians with disorder defined for cubic lattice clusters in 
three dimensions (3D) with many possible BC and inte- 
grate over all of them. We shall also consider the simpler 
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quantum crossover from the ballistic to localized phase in 
ID disordered systems. In this case we obtain a univer- 
sal Gaussian-like distribution at the crossover where the 
scaling parameter of the localization length £ over size L 
is fixed to unity. 

Quantum chaotic systems are well described by ran- 
dom matrix ensembles via Random Matrix Theory 
(RMT) HHIi]. The ievel-spacing distribution for small 
S involves the level-repulsion condition P(S) oc S 13 , 
where the so-called universality class index (3 — 1,2,4 
is determined from the remaining fundamental symme- 
tries of quantum chaotic systems. In the presence of time 
reversal invariance one obtains the orthogonal universal- 
ity class with [3=1. An applied magnetic field which 
breaks the time-reversal invariance leads to unitary uni- 
versality class with (3 — 2 while spin-orbit coupling with 
time-reversal to symplectic universality class with (3 = 4. 
The critical P(S) is also not independent of the remain- 
ing symmetries of a chaotic system either so that appro- 
priate expressions for the corresponding critical P(S) in 
different universality classes can be found [lflj . In the rest 
of the paper we shall rely on numerical tools in order to 
explore level statistics at criticality. Our aim is to answer 
the previously posed question whether the semi-Poisson 
distribution appears after averaging over all possible BC. 

We obtain spectra and their fluctuations in cubic 3D 
lattices with disorder. We have used standard diagonal- 
ization algorithms to compute eigenvalues for complex 
Hermitian sparse random matrices of size L 3 which in- 
clude flux. Many possible BC are considered by adding 
fluxes in all directions. More details of the numerical 
BC integration method can be found in where it was 
used for Hubbard models to describe interacting elec- 
trons. In the quantum coherent chaotic metallic regime 
the results for the eigenvalue level-spacing distribution 
P(S) are described by the Wigner surmise In the 
insulating integrable regime quantum coherence is lost 
since the localized states do not communicate with each 
other located in random positions, so that one obtains 
the random Poisson distribution. At criticality after in- 
tegration over all BC we obtain a critical P(S) which is 
closely described by semi-Poisson. However, we also find 
rather small deviations, particularly in the large S tail of 
the distribution. 

We consider a tight-binding model with site random- 
ness defined on a L x L x L finite lattice described by 
the Hamiltonian with potential energy and kinetic en- 
ergy terms 



H 



(cJjC m + cJ„C n ). 

(nm) 



(1) 



ri j cjj is a corresponding creation operator of an elec- 
tron and V n is the random site potential, a real inde- 
pendent random variable satisfying a box distribution 
P(Vn) = w> f° r — T" — < \ of zero mean and 
width W which measures the strength of disorder. We 
impose the generalized periodic boundary conditions on 
the wave functions ip(x, y, z) in all directions 



ip(x + L,y,z) = e %a *il>(x,y,z) 
^(x,y + L,z) = e ia ^ip(x,y,z) 
ip(x, y,z + L) = e lolz ip(x, y, z), 



(2) 



The sum is taken over all lattice sites n and all bonds 
(nm) where n, m denote nearest neighbor lattice sites, 
c„ is an annihilation operator of an electron on site 



which is equivalent to applying Aharonov-Bohm mag- 
netic flux a with components (a x ,a y ,a z ) varying in 
[0,2tt). 

For the zero disorder W = ballistic limit Eq. (1) 
has eigenvalues €j xt j j M (k x , k y , k z ) = — 2(cos k x +cos k y + 
cosfc z ) and plane wave eigenstates l >pj x ,j y ,j z (x,y, z) = 
L~ 3 / 2 exp(— i(k x x + k y y + k z z), with the positions of the 
A;-vectors determined by the BC of Eq. (2) via 

/■'- "' ■ jx = 0,..,L-l (3) 

_ a y + 2wjy 

Ky — — , Jy — U,...,±J 1 

a z + 2irj z 
k z = ^ , h = 0, ...,L- 1. 

For a given set of BC the allowed values of k x ,k y ,k z 
form a rigid grid of L x L x L points in the Brillouin 
zone. By varying the a x . ytZ s within [0, 27r) each Appoint 
in the Brillouin zone shifts to cover its own box so that 
the boxes of all points partition exactly the Brillouin zone 
When Oi x ,y,z = or 7r the added flux is equivalent 
to periodic, antiperiodic BC in the appropriate x, y, z 
direction, respectively. For critical disorder W c = 16.4 
0,0 Eq- (1) nas eigenvalues £j m j y ,j z and fractal eigen- 
states if>j tt j t j t (x, y,z). For any disorder W and size L the 
eigenvalues are obtained by numerical diagonalization. 

To characterize spectral fluctuations we compute the 
nearest level-spacing P(S) distribution function. First, 
we make the local mean level-spacing A(E) oc l/p(E) 
constant equal to one where p is the mean density of 
states, by "unfolding" the spectrum via local rescaling of 
the energy. For this purpose one needs the local average 
level-spacing A(E) which can be obtained over many lev- 
els around E. Alternatively, the ith "unfolded" energy 
level Si — Afav(Ei) can be obtained from the averaged 
integrated spectral density Af av (E), leading to average 
nearest level-spacing (S) = (£; — = 1. In the stud- 

ied energy regime around the band center the results did 
not change much with or without a proper unfolding since 
the mean density of states is almost constant 7]. 

In Fig. 1 we present the obtained level-spacing distri- 
bution P(S) in 3D for various sizes L with disorder W 
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FIG. 1: Metal-insulator transition. The level-spacing distri- 
bution function P(S) is shown in 3D disordered systems for 
the integrable, critcal and chaotic regimes for 10 random con- 
figurations of linear size L = 4, 10 by integrating over 1000 
possible boundary conditions, including periodic and antiperi- 
odic, corresponding to 1000 points in fc-space making up a 
total number of a few million eigenvalues. As the system 
size increases the data for the metallic regime with W — 10 
and the insulating regime with W = 30 approach the dashed 
lines on the right and left of the figure, which are the Wigner 
surmise P(S) = (w/2)S exp(— (tt/4)S 2 ) and the Poisson law 
P(S) = exp(— S), respectively. The filled dots of the com- 
puted scale-invariant distribution at critical disorder strength 
W c = 16.4 are compared with the semi-Poisson distribution 
P(S) = 45'exp(-2S') denoted by the solid line. 



by integrating over 1000 boundary conditions (BC). We 
have used the majority (about a third of the total) of en- 
ergy levels centered around the band center ignoring only 
those at band tails. For a given finite size L the results of 
disorder strength W corresponding to localization length 
£ lie in-between the Wigner surmise (L — ► oo asymptotic 
for extended states) and the Poisson distribution (L — > oo 
asymptotic for localized states) Q • Only for large enough 
L 3> £ the data eventually approach the Poisson curve 
which for W > W c indicates only localized states in 3D. 
For W < W c the data lie close to Wigner surmise which 
is also quickly approached for large size. Our results of 
Fig. 1 below the critical point (W — 10 < W c ) demon- 
strate a P(S) which approaches Wigner and above the 
critical point (W — 30 > W c ) a P(S) which approaches 
Poisson, respectively. At W — W c the obtained scale- 
invariant critical P(S) averaged over BC has Wigner-like 
small spacing behavior and Poisson-like large spacing be- 
havior. It is roughly described by the semi-Poisson curve 
although some differences become obvious, particu- 
larly in the tail of the distribution as seen in the semi-log 
plot of Fig. 2. These deviations also seem to violate the 
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FIG. 2: A linear-log plot of the critical P[S) shown in Fig. 
1 where prominent deviations from the semi-Poisson can be 
seen in the large-i? tail. 



scale-invariance of the critical curve, since they slightly 
increase by increasing L. However, in this part of the 
distribution P(S) is very low and it is much more likely 
that numerical accuracy of the data is lost. 

We have also examined the simpler ballistic to local- 
ization crossover. It occurs in ID disordered systems be- 
tween two integrable limits, one purely quantum where 
the kinetic energy term of H in Eq. (1) dominates 
and the other classical characterized the fluctuating ran- 
dom potential energy term of Eq. (1). In the ballistic 
regime one asymptotically reaches a simple ^-function 
P(S) while in the localized regime the Poisson distribu- 
tion. Our results taken at the band center E m are 
displayed in Fig. 3. Since the localization length is eas- 
ily computed with the anomalous perturbational result 
£ w 105.24H/- 2 at E = [HEI, di fferent from usual 
perturbation £ ks 96W~ 2 ^l — (E/2) 2 valid at most en- 
ergies E. We vary W for many random configurations 
which changes £ by considering appropriate sizes L to 
keep the scaling ratio £/L fixed. In the ballistic limit 
we find that a delta-function P(S) is approached for 
£/L ^> 1 while in the localized limit Poisson is reached 
for ijh < 1. For f/Z = 1/5 we obtain a P(S) closer to 
Poisson shown in Fig. 3(b). In Fig. 3(a) we display the 
obtained crossover distribution at the fixed ratio £/Z = 1 
for various £'s (varying W) and Z's. We observe that the 
scale-invariant crossover distribution is similar to a sim- 
ple Gaussian. Of course differences exist, such as the 
requirement of positive S which does not allow the exis- 
tence of Gaussian tails moving in the left-hand side. We 
emphasise that the computed P(S) measures fluctuations 
over the disorder ensemble precisely at the band center, 
rather than fluctuations in the energy domain studied in 
billiards. 
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FIG. 3: Ballistic to localized crossover. The level-spacing 
distribution function P(S) in ID disordered systems is shown 
to depend on the scaling parameter where £ is the lo- 

calization length and L the size, (a) In the crossover regime 
between ballistic and localized states having fixed £/L = 1 
from diagonalization of 100000, 500000 and 1000000 random 
matrixes of size L — 500, 1000, 5000 for energies close to E = 
with appropriate disorder values W. (b) In the localized 
regime having £/L = 1/5 from diagonalization of 500,000 ran- 
dom matrices of L = 500, 1000, 5000. In the ballistic regime 
£/L > 1 the P(S) (not shown) rapidly approaches a delta 
function centered around the mean 5 = 1. 

We have presented numerical results for the level- 
spacing distribution P(S) of non-interacting electrons in 
the presence of potential energy fluctuations at the An- 
derson metal-insulator transition. The critical behav- 
ior characterized by multifractal eigenstates leads to a 
generic sensitivity to BC In the present work im- 
posed Ahai'onov-Bohm fluxes are equivalent to consid- 
ering many BC. For example, a continuous set of BC 
in ID leads to the equivalent problem of studying the 
non-overlapping mini-band widths [14j. The AB-flux 
controlled crossover from orthogonal to unitary studied 
in has also revealed sensitivity to the imposed AB 
fluxes. Although such sampling over the Brillouin zone 
via the Aharonov-Bohm fluxes breaks time-reversal in- 
variance, by changing the universality class from orthog- 
onal to unitary, the obtained results remain valid for the 
orthogonal case. 

In summary, the P(S) distribution for disordered sys- 
tems is computed to test the semi-Poisson form at crit- 
icality. The obtained results complement previous stud- 
ies for understanding the metal-insulator transition. Our 
main finding is that the corresponding critical statistics 
is reasonably well described by the semi-Poisson function 
after an average over many BC is taken. Small but rather 
persistent deviations occur mostly in the difficult to reach 
large S tail, possibly due to numerical artifacts arising 
from limited accuracy. The spectral fluctuations stud- 



ied here could help to understand various properties of 
current-currying states in mesoscopic systems and might 
also suggest a critical RMT. A disorder system is seen 
to progress from ballistic extended (periodic integrable) 
to diffusive (chaotic) and finally to localized (random in- 
tegrable) regimes. The ballistic regime is dominated by 
quantum correlations and the localized regime by random 
potential energy fluctuations, respectively. The interme- 
diate quantum chaotic regime is reached either by shape 
peculiarities in billiards or for moderate disordered po- 
tential in solids. We add the remark that the Gaussian- 
like P(S) observed for the addition s pec tra of the con- 
ductance peaks in Coulomb blockade (lfij . if it could be 
interpreted within a non-interacting electron framework, 
would suggest an analogy to the obtained crossover be- 
tween two integrable regimes with and without disorder. 

I like to thank H. Schomerus for a critical reading of 
the manuscript. 
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